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Abstract. The Milky Way has been estabished to emit gamma rays. These gamma rays are presumably dominated by decays 
of neutral pions, although inverse Compton scatterings and bremsstrahlung also contribute. It is plausible that other galaxies 
can be diffuse sources of gamma rays in a similar manner. Starburst galaxies are particularly interesting to study as they are 
expected to have much higher cosmic-ray fluxes and interstellar matter densities. The neutral pions are created in cosmic-ray 
interactions with interstellar matter. Presented here is an overview of the recent work by Karlsson and co-workers on proton- 
proton interactions and the resulting secondary particle inclusive cross sections and angular distributions. This model can be 
used to calculated the ;r" component of the gamma-ray yield and spectrum from a starburst galaxy. The yield is expected to 
increase significantly (30% to 50%) and the spectrum to be harder than the incident proton spectrum. 
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INTRODUCTION 

The Milky Way Galaxy is a known source of gamma 
rays in MeV to GeV energies. This emission is partic- 
ularly strong in the Galactic disk and traces closely the 
distribution of gas in the Galaxy. It is presumably dom- 
inated by gamma rays from decaying neutral pions, cre- 
ated in interactions between cosmic rays and interstellar 
matter (see e.g. The gamma-ray flux and spectral 
shape measured by the Energetic Gamma-Ray Experi- 
ment Telescope (EGRET) |2] is considered the key at- 
testation of this interpretation. Contributions also come 
from inverse Compton scatterings and bremsstrahlung 
due to electrons. 

It is plausible that other galaxies with properties sim- 
ilar to the Milky Way also emit gamma rays, but the ex- 
tremely large distances involved and low expected lumi- 
nosities make detections very unlikely. To date, the Large 
Magellanic Cloud is the only external galaxy detected as 
a diffuse source of gamma rays ll]. 

In a starburst galaxy there are regions with star- 
formation rates several orders of magnitude higher than 
that of normal galaxies. With an enhanced star-formation 
rate follows an enhanced rate of supernovae explosions 
and the cosmic -ray flux is expected to be higher (un- 
der the assumption of cosmic-ray acceleration in super- 
nova remnants). Considering a galaxy with much higher 
cosmic-ray flux and a denser interstellar medium, com- 
pared to the Milky Way, the gamma-ray luminosity may 
be high enough to allow for a detection, even in TeV 
gamma rays, despite the distances involved. 

Presented here is the parameterization of gamma-ray 
inclusive cross sections |^ and angular distributions tSJ 



produced by proton-proton interactions in astronomical 
environments. The parameterization is derived from sim- 
ulations using the Pythia Monte Carlo event generator 
The parametric model can be used to model the hadronic 
component of the gamma-ray emission from a typical 
starburst galaxy and allows accurate calculations of the 
expected gamma-ray flux and spectral shape. The fluctu- 
ations in the histograms are due to low event statistics. 

HADRONIC PRODUCTION OF GAMMA 
RAYS 

Pions are produced in hadronic interactions between cos- 
mic rays and interstellar matter. The major contribution 
comes from proton-proton (p-p) interactions, but there 
are also contributions from heavier nuclei. The neutral 
pion, 71^, decays predominantly into two gamma rays and 
the charged pions decay into muons, which subsequently 
decay into electrons/positrons, and neutrinos. 

The work by Kamae et alj_||6|] showed that previous 
models (see references in fsf]) of pion production left 
out two important features of the inelastic p-p interac- 
tion; the diffraction dissociation process and the Feyn- 
man scaling violation. It was also noted that the inelastic 
p-p cross section was assumed to be constant, about 24 
mb, for Tp 10 GeV, in contradiction to recent exper- 
imental data, where a logarithmic increase with proton 
energy is evident This is shown in Figure [T] together 
with the component cross sections used in this work. 

The model by Kamae et al. |6] included only the 
non-diffractive interaction and the diffraction dissocia- 
tion process. This model is not accurate near the pion 
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FIGURE 1. Experimental p-p (p-p) cross sections, total 
(squares) and elastic (triangles), as functions of proton labo- 
ratory momentum, and their model fits; total (thin solid) and 
elastic (dot-dot dashed). The other lines are: total inelastic 
(thick solid), non-diffractive (dashed), diffractive (dot-dashed), 
A(1232) (dotted), and res(1600) (thin dashed). 

production threshold (see Figure 5 in Kamae et al. |0]). 
To correct this and improve accuracy the contribution 
of two baryon resonance excitation states, the physical 
A(1232) resonance and res(1600) representing several 
resonances around 1600 MeV/c^, were included. Figure 
|2] shows the improved total inclusive %^ cross section. 
Further necessary adjustments to model A required by 
the baryon resonances are described in section 3 of Ka- 
mae et al. 

The interaction model presented here does not include 
a-p, He, and a -He due to paucity of experimental 
data. To an approximation, the a-particle and the He 
nucleus can be taken as four individual nucleons. The 
error is expected to be less than 10% for high-energy 
gamma rays. 



Monte Carlo Simulations 

The parameterizations of inclusive cross sections and 
angular distributions are derived from Monte Carlo sim- 
ulations of the adjusted p-p interaction model described 
in Kamae et al. f?]. Events were generated for each of 
the four components at discrete proton kinetic energies 
(0.488 GeV < 7), < 512 TeV) from a geometrical series 

Tp = 1000 X 2('-22)/2 GeV, / = 0, . . . ,40. (1) 

Each proton kinetic energy, Tp, represents a bin covering 
2-0.25 to The addition of the resonances to the 
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FIGURE 2. Experimental (squares) and Monte Carlo simu- 
lated (solid) total inclusice ifi cross section as functions of pro- 
ton laboratory momentum. The other lines are for the four com- 
ponents of the p-p interaction model: non-diffractive (dashed), 
diffractive process (dot-dashed), A(1232) (double-dot dashed), 
and res(1600) (dotted). 



model required an increased sampling frequency near the 
pion production threshold and events were also generated 
for Tp — 0.58 and 0.82 GeV. Events were not generated 
for proton energies where the component cross section is 
very small or zero. 

The results of these Monte Carlo simulations have 
been verified to agree with experimental data available 
forpions. The total inclusive tt" cross section agrees with 
the experimental one (see Figure ^ and the simulations 
were ensured to reproduce the distributions of pion ki- 
netic energy in the p-p center-of-mass (CM) system in 
the resonance region (see Figures 3,4 and 5 in [4]). The 
angular distribution was verified using experimental data 
on the Lorentz invariant cross section, EcPa /dp^, and 
the average transverse momentum of pions (see Figures 
2 and 3 in lH). 



Parameterizations 

The inclusive cross sections, A(7(£',rp)/Alog£', of 
all stable secondary particles, gamma rays, electrons, 
positrons, electron and muon neutrinos (and their anti- 
particles) have been parameterized. For each component, 
the mono-energetic cross section Act (/s) /A log/: were 
fitted to a common algebraic form F{x) with x = \ogE 
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FIGURE 3. 



Gamma-ray spectram produced by protons with 
a power-law spectrum of index 2.0 and cutoff at Ty, = 512 
TeV. The solid line was calculated using the parameterization 
described in the text and the histogram was derived from the 
Monte Carlo simulated events. The dashed line is an eyeball fit 
to a power law with index 1.95. The absolute normalization is 
relative to the density of target protons. 




FIGURE 4. Gamma-ray spectra from a pencil beam of pro- 
tons, power law with index 2.0, observed at three different an- 
gles: head on (solid), 0.5° (dashed) and 2° (dotted). The thin 
line is for integration over the entire phase space. Histograms 
are the corresponding spectra from Monte Carlo simulations. 
The absolute normalization is relative to the density and distri- 
bution of target protons. Fluctuations in the histograms are due 
to low event statistics. 



and E in GeV. Each fit gives a set of parameters which 
were then fitted as functions of Tp to give the final param- 
eterization of the inclusive cross section. The algebraic 
forms and the fits are described in Kamae et al. jj] . 

The gamma-ray spectrum, E^dF /dE, produced by a 
power-law distribution, index 2.0, of protons was calcu- 
lated using the parameterization and compared with the 
corresponding Monte Carlo spectrum in Figure |3] The 
parameterization is accurate to within 10% except near 
the upper kinematic limit where the error is up to 20%. 

The angular distribution, i.e. the distribution over the 
polar angle 9, is readily derived from the pt distribution. 
The latter is expected to follow a simple exponential 
form. Thus, for gamma rays, the pt distribution was fitted 
to an exponential for each Tp and also the gamma-ray 
energy E. Again, each fit gives a set of parameters which 
were fitted as functions of Tp and E to give the final 
parameterization of the distribution. The details of this 
procedure are described in Karlsson and Kamae |5]. 

The gamma-ray spectrum due to a beam of protons 
along the z-axis with no spatial extension in the x-y 
plane, i.e. a pencil beam, has been calculated using the 
parameterization of p, distributions. The beam is viewed 
from three different angles, head on, 0.5° and 2°, giving 
three different spectra. These is shown in Figure |4] to- 
gether with the total inclusive spectrum integrsted over 
the entire phase space. Each spectrum is integrated over 
the annular portion (9,9 + A9) with A9 — 2'. The nor- 
malization is relative to the density and distribution of 
target material. The histogram fluctuations are due to low 
event statistics. 



STARBURST GALAXIES 

A starburst galaxy is a galaxy with regions on highly en- 
hanced star formation. Compared with a normal galaxy, 
the star-formation rate can be orders of magnitude higher 
in a starburst region. This also implies a higher rate of 
supernovae explosions and if the assumption that super- 
novae remnants accelerate the bulk of the cosmic rays 
in a galaxy is correct, then the total cosmic-ray flux in a 
starburst region is expected to be much higher than what 
is observed in the Milky Way. A starburst region also 
has a denser interstellar medium. This make starbursts 
the more likely candidates for detection of a non-active 
galaxy, even in TeV gamma rays. 

Several predictions of expected gamma-ray fluxes 
from starburst galaxies have been made. Paglione et al. 
Q modeled NGC 253 and matched their results to data 
obtained with EGRET. The EGRET data gives a 2(7 up- 
per limit on the flux above 100 MeV of 8 x 10"** photons 
cm"^ s"^ Based on their calculations of the gamma-ray 
spectrum (see Figure 6 in |8]), the Fermi Large Area 
Telescope (Fermi-LAT) should be able to detect NGC 
253 based on a one-year observation. Current generation 
ground based telescope will most likely have difficulties 
detecting this source. A detection was claimed by CAN- 
GAROO |9] but the spectrum is uncertain and H.E.S.S. 
collaboration has not confirmed the detection. 

Perhaps an even better detection candidate is M82 
in the northern hemisphere. The integrated flux above 
100 GeV was estimated to be about 2 x 10"'^ photons 
cm"^ s"' ifioll . Based on this model, it is a possible 



candidate for detection by the ground-based telescopes 10. M. Persic, Y. Raphaeli, and Y. Arieli, AcfeA 486, 143-149 
MAGIC and VERITAS and it is a good candidate for the (2008). 
MAGIC II telescope. They also calculated the integrated 
flux above 100 MeV to be about xlO^^ photons cm^^ 
s^^ which is detectable by the Fermi-LAT over a one- 
year observation. 

The author has noted that most calculations of the pion 
component of the gamma-ray spectrum is based on the 
outdated p-p interaction models discussed earlier in this 
paper In particular, they assume a constant inelastic p-p 
cross section. As was shown by Kamae et al. inclu- 
sion of a logarithmically rising cross section significantly 
increases the gamma-ray yield at TeV energies. 

CONCLUSIONS AND FUTURE WORK 

The p-p interaction model presented here facilitate com- 
putation of gamma-ray spectra in both cases of both 
isotropic and anisotropic proton distributions. The for- 
mulae incorporate all important known features of the in- 
elastic p-p interaction up to about 500 TeV. Its inclusion 
in starburst models is expected to increase the gamma- 
ray yield by at least 25%. 

The author aims to further refine the p-p interaction 
model and apply it to detailed modeling and analysis of 
several starburst galaxies, including prominent detection 
candidates M82 and NGC 253. A more rigorous study of 
several other starbursts will be earned out. 
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